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Abstract. With the recent advances in radioactive ion beam technology, Coulomb excitation at safe energies 
becomes an important experimental tool in nuclear-structure physics. The usefulness of the technique to 
extract key information on the electromagnetic properties of nuclei has been demonstrated since the 1960’s 
with stable beam and target combinations. New challenges present themselves when studying exotic nuclei 
with this technique, including dealing with low statistics or number of data points, absolute and relative 
normalisation of the measured cross sections and a lack of complementary experimental data, such as 
excited-state lifetimes and branching ratios. This paper addresses some of these common issues and presents 
analysis techniques to extract transition strengths and quadrupole moments utilising the least-squares fit 
code, GOSIA. 

PACS. 25.70.De Coulomb excitation - 21.10.Ky Electromagnetic moments - 29.38.Gj Reaccelerated 
radioactive beams - 29.85.Fj Data analysis 


1 Introduction 

Recent advances in radioactive ion beam (RIB) technol¬ 
ogy, in particular the increasing range of species and post¬ 
acceleration energies available from ISOL facilities such as 
REX-ISOLDE at CERN, SPIRAL at GANIL and ISAC 
at TRIUMF, has led to a resurgence of the use of nuclear 
reactions to study the structure of nuclei [1]. Specifically, 
Coulomb excitation at safe energies with RIBs is now giv¬ 
ing us a wide range of information on the electromagnetic 
properties of exotic nuclei, leading to knowledge of the 
nuclear shape or, more precisely, nuclear charge distribu¬ 
tion [2]. 

“Safe” Coulomb excitation is defined as the process 
of inelastic scattering of nuclei via the electromagnetic 
force such that the energy in the centre-of-mass frame en¬ 
sures negligible contribution to the reaction process from 
the strong force. This is fulfilled by maintaining a min¬ 
imum distance of 5 fm between the nuclear surfaces, of¬ 
ten called Cline’s “safe energy” criterion [3]. Exploiting 
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the well-understood electromagnetic interaction allows a 
nuclear-model-independent interpretation of the observed 
data. With the use of light-ions as probes, the excita¬ 
tion modes are often limited to single transitions from 
the ground state. This data can be interpreted in terms 
of a semi-classical description using first-order perturba¬ 
tion theory. However, the use of high -Z probes has meant 
that multiple-step excitation is now common, and a large 
number of states can be accessed from ground or isomeric 
states. The technique of data analysis based on coupled- 
channel calculations with the GOSIA code [4] has allowed 
for the determination of large, and in some cases com¬ 
plete, sets of low-lying E 2 and E 3 matrix elements in 
multi-step Coulomb-excitation experiments, including di¬ 
agonal matrix elements related to the static electromag¬ 
netic moments. Due to this completeness of measurement, 
low-energy Coulomb excitation with heavy ions (or high -Z 
targets) is an extremely sensitive probe of collective nu¬ 
clear structure. Used in conjunction with complementary 
spectroscopic data, such as excited-state lifetimes, y-ray 
and conversion electron branching ratios, multipole mix¬ 
ing ratios, electric and magnetic moments, mean-square 
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charge radii etc., a pure experimental understanding of 
low-lying collective modes and shapes can be achieved. 

When studying exotic nuclei with this technique how¬ 
ever, new challenges emerge. These include dealing with 
low statistics and a lack of complementary experimental 
data such as excited-state lifetimes and branching ratios. 
For many short-lived nuclei, especially on the neutron-rich 
side, precise information on the lifetimes of excited states 
is not known and thus another solution for the normali¬ 
sation of the measured Coulomb-excitation cross sections 
needs to be applied. In general two options are possible: 
either normalisation to the excitation of target nuclei with 
known electromagnetic matrix elements or to the number 
of elastically scattered beam particles. 

This paper attempts to address some of the common 
problems and solutions encountered with the extraction 
of electromagnetic matrix elements from RIB Coulomb- 
excitation experiments in general, with examples taken 
from studies performed at REX-ISOLDE and GANIL. 
Here, the GOSIA code (see Section 2) is most commonly 
used for this purpose. Firstly though, the observables from 
such experiments must be clearly defined; this is done in 
Section 3. Methods utilising the GOSIA code for the analy¬ 
sis are presented in Section 4 and a summary and outlook 
is given in Section 5. 

2 The gosia code 

Experiments performed in the 1950’s utilising light-ion 
beams as a means of exciting target nuclei were relatively 
simple to interpret using first- and second-order perturba¬ 
tion theory. Later, heavy-ion beam experiments populated 
many excited states via multiple-step Coulomb excitation. 
Early versions of computer codes designed to handle the 
analysis of these data, most notably that of Winther and 
de Boer [5], employed the semi-classical theory of multiple 
Coulomb excitation developed by Alder and Winther [6]. 
This code allowed quantitative calculations of excitation 
amplitudes for the first time, using a set of reduced electro¬ 
magnetic matrix elements as input. With this philosophy, 
the GOSIA code [4] was designed in 1980 to achieve an ex¬ 
traction of the electromagnetic matrix elements from a set 
of Coulomb-excitation data by performing a fitting routine 
using these matrix elements as parameters. Both excita¬ 
tion and the consequent y-ray de-excitation, governed by 
the very same set of matrix elements, are calculated within 
the code, allowing for a direct comparison to experimental 
data [4]. The description of y-ray de-excitation in GOSIA is 
based on the CEGRY code [4, 7] and takes into account the 
angular correlations, deeorientation effect, recoil effects, 
Jacobian to a common reference frame, and integration 
over detector geometry. 

The first successful application of GOSIA was to prove 
that the set of matrix elements obtained for 110 Pd [8] con¬ 
stituted a unique solution, which has been later confirmed 
by the results of a recoil-distance lifetime measurement [9] . 
The GOSIA code was further validated and tested in the 
analysis of an extensive data set for several W-Os-Pt iso¬ 
topes [10-13]. In the following years, it was used to study 


shape evolution and coexistence in many regions of the 
nuclear chart, including transitional nuclei [14-16], rare 
earths [17] and actinides [18], as well as exotic octupole 
shapes [19]. 

3 Observables in Coulomb-excitation 
experiments 

The direct observables in Coulomb-excitation experiments 
are usually the y-ray intensities corresponding to the scat¬ 
tering of the projectile particle defined by the observation 
of at least one of the collision partners in a given angu¬ 
lar and energy range. In contrast, the deduced matrix el¬ 
ements are not direct observables and usually occur as 
strongly correlated parameters in a fit of the y-ray inten¬ 
sity data. In order to relate these gamma-ray intensities to 
the excitation cross sections of the populated states, which 
can be calculated for a given set of scattering and nuclear 
parameters, normalisation factors need to be introduced 
as described in Section 4. 

Data sets introduced to GOSIA are most often described 
in terms of “experiments”. These may be defined by dif¬ 
ferent combinations of beam and target, beam energy and 
scattering-angle range. With the use of segmented parti¬ 
cle detectors, such as the Double-Sided Silicon Strip De¬ 
tectors (DSSSD) or Parallel Plate Avalanche Counters 
(PPAC), subdivision of the data can be made in terms 
of scattering angle, gaining sensitivity to second-order ef¬ 
fects such as the spectroscopic quadrupole moment, Q s . 
This is demonstrated in Figure 1 where the reorientation 
effect [20, 21] leads to an increasing deviation in the cross- 
section at large scattering angles for the assumption of dif¬ 
ferent quadrupole moments. This can be further increased 
by the use of different targets to disentangle contributions 
from single- and multiple-step excitation processes. 

It should be noted that in contrast to other spectro¬ 
scopic methods Coulomb excitation is not only sensitive to 
magnitudes of the electromagnetic matrix elements, but 
also to their relative signs having a direct influence on 
excitation probabilities. As an example one can consider 
a state A that can be populated in one-step E 2 excita¬ 
tion from the ground state or in a two-step E 2 excitation 
process via an intermediate state B , as depicted in Fig¬ 
ure 2(a). For each of the two possible excitation paths the 
contribution to the total excitation amplitude is propor¬ 
tional to the relevant matrix elements: (A||E2||g.s.) for 
the direct excitation and the product of (A||E2||11) and 
(R||l?2||g.s.) for the two-step process. The excitation 
probability is proportional to the square of the sum of 
excitation amplitudes and therefore it contains not only 
quadratic terms ((A||E2||g.s.) 2 , related to B(E2; A —> g.s.), 
and (A||E , 2||R) 2 (l?||E2||g.s.) 2 ) but also interference 
terms between possible excitation paths, such as 
(A||E2||g.s.)(A||£’2||13)(13||E2||g.s.). The signs of these in¬ 
terference terms depend on the relative signs of the matrix 
elements. This is illustrated by the example of 110 Ru [level 
scheme shown in Figure 2(b)] on 208 Pb presented in Fig. 3 
where for large scattering angles the population of the 2^ 


M. Zieliriska et al.: Analysis methods using the GOSIA code 


3 



Centre-of-mass scattering angle [deg] 

Fig. 1: Coulomb-excitation probabilities (dashed lines) 
and cross sections (solid lines; product of the Ruther¬ 
ford cross section and Coulomb excitation probability) for 
populating the 2f state in 184 Hg incident on a 120 Sn at 
2.8 MeV/u under the different assumptions for the spec¬ 
troscopic quadrupole moment, Q s . The oblate (black) as¬ 
sumption is that of Q s = 1.15 eb, extracted from the mea¬ 
sured B(E2;2^ —> Of) [22, 23] and the rigid-rotor model 
(K = 0), while the prolate assumption (red) has the same 
magnitude, but a negative sign for Q s . The spherical as¬ 
sumption (Q s = 0 ) is shown in blue. 



Fig. 3: Relative population of excited states in 110 Ru 
Coulomb excited on 208 Pb at 430 MeV beam energy, cal¬ 
culated for two different signs of ( 2 f ||l? 2 || 2 f): negative 
(solid lines) and positive (dotted lines), while all other 
matrix elements remain the same. The Of ground state is 
not shown, but dominates the remainder of the population 
at all angles. 

4 Coulomb-excitation data analysis 

4.1 Normalisation of measured cross sections 
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Fig. 2: (a) Schematic level scheme showing the various 
excitation paths that lead to sensitivity to the relative 
signs of matrix elements, see text for details, (b) Low-lying 
states in 110 Ru included in calculations shown in Fig. 3. 


In order to extract nuclear-structure parameters (matrix 
elements) from Coulomb-excitation data, the measured 7 - 
ray intensities have to be converted to absolute excitation 
cross sections. Possible complications arise from the fact 
that the efficiency of the particle detection set-up, dead¬ 
time, beam intensity etc. are not always known with good 
precision. To deal with this, GOSIA uses normalisation con¬ 
stants, which relate the calculated and experimental in¬ 
tensities. These can be fitted or given by the user, as de¬ 
scribed in the following sections. In the most general form, 
the normalisation constant used in GOSIA is a product of 
the Rutherford cross section, the time integrated beam 
current, the absolute efficiency of particle and 7 -ray de¬ 
tection and the particle solid angle factor. If the statistics 
are not sufficient to make use of particle- 7 -ray angular cor¬ 
relations (which is usually the case for radioactive beam 
studies), 7 -ray spectra from individual detectors may be 
summed together, reducing the number of necessary nor¬ 
malisation constants to one per experiment. In such cases, 
the relative 7 -ray detection efficiency as a function of en¬ 
ergy has to be provided for each detector. 

The normalisation constant, C, for a given experiment 
is fitted to all measured 7 -ray intensities I e observed in 
an experiment by minimising the expression: 


state depends very strongly on the sign of ( 2 (j t_ ||£ 12 || 2 ^") 
with respect to those of ( 2 ] I "||.E 2 || 0 ] , ~) and ( 2 ^"||i? 2 || 0 ] 1 "). 
This effect can be strong enough to be visible even in low- 
statistics RIB measurements and thus for several exotic 
nuclei complete sets of matrix elements including their 
relative signs have been determined [23-26]. 


J2(C^ - i?) 2 /af (1) 

i 

where If denotes the calculated 7 -ray intensity for the i-tli 
observed transition integrated over beam energy and scat¬ 
tering angle, If its measured intensity and cq its experi¬ 
mental uncertainty. Moreover, it is possible to introduce 
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relative normalisation constants C m that link data sets re¬ 
sulting from the subdivision of data collected during one 
physical run into m slices of scattering angle. If for each of 
m coupled experiments a relative normalisation constant 
Cm is defined in GOSIA, during the minimisation of the y 2 
function the following expression is minimised and only 
one global normalisation constant C g i 0 b a i is fitted: 

Cnlt-Itf/a? ( 2 ) 

m i 

It should be noted that the C m factors can be arbitrarily 
rescaled, as the scaling factor can be always incorporated 
in Cgiobai- The normalisation constants are fitted in GOSIA 
at the same time as the matrix elements, during the min¬ 
imisation of the x 2 function described in Section 4.2. 

The products C m If that are compared to experimen¬ 
tal 7 -ray intensities depend obviously both on the matrix 
elements and on the normalisation constants. Especially 
in the case of one-step excitation, one can easily com¬ 
pensate a modification of the relevant matrix elements 
by adjusting the normalisation constant. Therefore, in or¬ 
der to obtain a reliable set of matrix elements, additional 
constraints on either the matrix elements or the normal¬ 
isation constants have to be provided. The possible tech¬ 
niques, depending on the specifics of the experiment, are 
presented in the following sections. 

4.1.1 Elastic scattering 

Historically, the simplest and most direct method of nor¬ 
malising Coulomb-excitation cross sections is to use the 
measured elastic-scattering (Rutherford) cross section. 
This requires precise knowledge of the scattering angu¬ 
lar range, efficiency of the particle detection system and 
well understood dead time if one is to obtain the inte¬ 
grated beam current. Since the Rutherford cross section 
is very sensitive to scattering angle at low centre-of-mass 
angles, uncertainties related to geometry are minimised 
for backscattering as demonstrated in Ref. [27, 28]. For 
inverse kinematics reactions, the backscattered projectiles 
are forward focused in the laboratory frame of reference 
and have low energy. The corresponding recoils however, 
can be utilised where clean kinematic separation of these 
events can be made. In RIB experiments, where beam in¬ 
tensities are low, the highest excitation probability is de¬ 
sired and as such, high -Z targets are usually used. Un¬ 
certainties are introduced because of events from different 
scattering angles that can be misinterpreted. This is par¬ 
ticularly true for many experiments utilising silicon strip 
detectors at forward angles, such as those at GANIL [24] 
and with Miniball at REX-ISOLDE [29]. 

In cases where absolute Rutherford cross sections are 
not reliable (for example when downscaling is applied to 
single particle events and thus the dead time is different 
as compared to particle-gamma events), particle singles 
events may still be used to calculate the relative normal¬ 
isation constants, C m . Commonly, this applies to experi¬ 
ments where data is taken in the same run but is divided 


into different angular cuts. For this, one needs knowledge 
of the total number of scattered particles in each angular 
range, N m , i.e. without a coincidence condition on 7 rays 
or a second particle. The two are related by the following 
expression: 

= a/I . 0) 

where {A9 m , represents the solid angle subtended in 
the experiment. Again, C m may be arbitrarily rescaled due 
to the remaining normalisation fitted by GOSIA, Cgiobai, 
but the ratios of each coupled C m remains the same. 

4.1.2 Excited-state lifetimes or B[E 2) values 

When multiple states are excited, with single- or multiple- 
step Coulomb excitation, one or more B{E2) values con¬ 
necting the ground-state and an excited state can be used 
to fit the normalisation constants for each experiment in 
GOSIA, C m . For this, one must also observe the corre¬ 
sponding population of such a state with good precision, 
which means that the 7 -ray intensity and efficiency, along 
with the branching ratio, must be known to good preci¬ 
sion. This is usually the simplest and preferred method in 
these cases as everything is fitted by the code and there 
are no additional calculations required by the user. 

In even-even nuclei, the normalisation is usually ful¬ 
filled by an independent measurement of the 2 j t “-state life¬ 
time, r(2f). Two examples of this technique with RIBs, 
are the cases of 74 ’ 76 Kr [24] and 182 - 188 Hg [23, 25, 30], 
where multiple lifetimes of yrast states were known in the 
literature and even re-measured [22, 31, 32] to provide 
the required precision. For odd-mass or odd-odd systems, 
multipole mixing ratios also become important since the 
strongest-observed 7 ray is often a mixed E2/M1 tran¬ 
sition (see also Section 4.4.3). Furthermore, low-energy 
transitions in heavy nuclei can also be strongly converted, 
meaning that the strongest excitation path may not result 
in an intense 7 -ray decay. In these cases, it is usually pos¬ 
sible to normalise to the next higher-lying transition since 
the low-energy of the first-excited state also means that 
the probability of two-step excitation approaches that of 
the single-step excitation, as was done in the analysis of 
224 Ra [33]. 

4.1.3 Target excitation 

The electromagnetic interaction between the collision part¬ 
ners causes excitation of either the projectile or target nu¬ 
cleus. The observed excitation of target nuclei can usually 
be described with high precision using literature values of 
relevant matrix elements and used to normalise the exci¬ 
tation cross sections measured for beam nuclei. The ob¬ 
served number of 7 rays in the transition de-exciting an 
excited state in the target nucleus, can be described in the 
following equation: 

N t = L ■ pd ^ A . fr t e 7 (E’t)ep art <7t (4) 

A-t 
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where cr t is the integrated cross-section of exciting the 
given state in the target, b t is the total 7-ray branching 
ratio for the transition, e^(E t ) is the absolute efficiency 
of detecting a 7 ray of energy E t , e part is the efficiency 
of detecting a particle in the angular range defined by 
the integration limits of the cross-section, pd is the thick¬ 
ness of the target in mg/cm 2 , N a is Avogadro’s number, 
A t is the mass number of the target and L is the time- 
integrated luminosity of the beam. A similar equation can 
be constructed for the number of 7 rays in the transition 
de-exciting an excited state in the projectile, assuming the 
same angular range for particle detection: 


Np — L • • bp€ry (Ep') ep ar t Op 

A t 


( 5 ) 


Taking a ratio of Equations 5 and 4 removes both the 
intrinsic particle detection efficiency and luminosity: 

Np _ b p e-y{Ep)ap . . 

N t ~ b t e 1 (E t )a t 1 j 


meaning that one can solve Equation 6 for a p and there is 
no requirement to have knowledge of the integrated beam 
current. This is the principle of GOSIA2. 

When dealing with RIBs, pure beams are often not 
achievable and the target is also excited by beam contam¬ 
inants. If the beam composition is monitored during the 
experiment, this can be dealt with rather simply with the 
following correction to the experimental 7-ray intensities 
from the target [ 34 ]: 


F = 


1 


1 + 



cr t(Z c ,A) \ ’ 
<r t (Z x ,A) ) 


( 7 ) 


where a t (Z,A) is the cross section of the target, excited 
by a beam with proton number Z and mass A. For every 
contaminant, c, with Z = Z c , the ratio to the component 
of interest with Z = Zx, can be expressed as r c = I c /Ix , 
where I Cj x is the intensity of the respective components 
in the beam. 

There also exists the possibility of impurities in the tar¬ 
get. In this case the experimental intensities measured for 
the beam must be corrected to account for the scattering 
on target impurities. For this, knowledge of the isotopic 
purity is required. This can be either from the target man¬ 
ufacturer or the observed excitation ratios, deduced from 
7-ray intensities, if available. Assuming only two compo¬ 
nents, a correction factor, Iq, can be calculated for each 
excited state, i [26]: 


Ft = 



*i{Z',A') \ 
&i(Z,A) ) 1 


( 8 ) 


where <Ji(Z, A) and ai{Z' , A') are the excitation cross sec¬ 
tions of a given state in the projectile on the main target 
species and contaminant, respectively. These can be cal¬ 
culated by GOSIA, obtaining the ratio for each transition 
given a set of starting matrix elements. The isotopic pu¬ 
rity, P , is expressed by 


P = 


N a 

N a P 


( 9 ) 


where Na,a< are the numbers of atoms of mass A, A'. By 
taking the ratio of the cross sections with different masses, 
at the same laboratory angles, the differences in Ruther¬ 
ford cross section and the centre-of-mass-dependent ex¬ 
citation probabilities are accounted for. However, F z re¬ 
mains an estimation since the excitation probability of 
each state will depend in a complex manner on the elec¬ 
tromagnetic matrix elements. A systematic error must be 
retrospectively estimated due to this assumption by re¬ 
calculating Fi with the final set of matrix elements. Dif¬ 
ferences between the original and final estimations of F t 
are usually small if P is large. In the case of 196 Po on 
94 ( 95 )Mo(P = 95(2)%) [26], the maximum systematic er¬ 
ror in Fi was calculated to be 0.6%, which is much smaller 
than the statistical uncertainty. 


4.2 x 2 square minimisation in GOSIA 

The set of electromagnetic matrix elements is extracted 
by performing the minimisation of the \ 2 function. The 
total x 2 function is built of measured 7-ray intensities and 
other known spectroscopic data, and those calculated from 
a set of matrix elements between all relevant states. The 
calculated 7-ray intensities are corrected for effects such 
as: internal conversion of electromagnetic transitions, the 
energy-dependent efficiency of the 7-ray detectors and the 
angular distribution of the emitted radiation. A proper re¬ 
production of the experimental 7-ray intensities requires 
integration over the scattering angular ranges, defined by 
the particle detection set-up, and over the range of inci¬ 
dent projectile energies resulting from the energy loss in 
a target. The convergence of the % 2 fit can be improved 
by using known spectroscopic data, e.g. 7-ray branching 
ratios, multipole mixing ratios or lifetimes. 

The x 2 function consists of three components resulting 
from various subsets of data: 

x 2 = S y + S l + S d . (10) 

The first contribution, S y , comes from the comparison 
of 7-ray intensities observed in the experiment, If., and 
those calculated with the fitted matrix elements, If, and 
is defined as: 

s y = Y J ^Y,^^io l k-n)\ (ii) 

ij k(ij) k 

The summations extend over all defined experiments, i, 7- 
ray detectors, j, and the detector- and experiment- depen¬ 
dent number of observed transitions indicated by k. The 
coefficients Cij are normalisation constants connecting ex¬ 
perimental and calculated intensities. These are equivalent 
to C m described in Section 4 . 1 , but the summation now 
extends independently over the number of independent 
7-ray detectors as well as experiments or sub-divisions. 
These are defined individually for each experiment and 
detector combination and fitted on the same basis as the 
matrix elements. The weights, Wij, ascribed to the various 
subsets of data defined by different experiments and 7-ray 
detectors, can be set independently by user. 
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The second contribution, Si, is related with the user- 
defined “observation limit” and is defined as follows: 



An experiment and detector dependent upper limit of 7 - 
ray intensities, u(i,j ), is expressed as a fraction of the 
normalising transition specified by the user (usually it is 
the strongest observed transition, i.e., 2 + —>• 0 ^ for even- 
even nuclei). If the calculated intensity of any unobserved 
7 -ray transition, divided by the intensity of the normal¬ 
ising transition, /“(i, j), exceeds this upper limit then it 
is included in the calculation of the least squares fit. The 
summation extends over the calculated 7 -ray transitions 
in each experiment and detector combination not defined 
as experimentally observed, provided that the upper limit 
has been exceeded. A proper set of upper limits prevent 
finding unphysical solutions yielding the production of 7 - 
ray transitions not observed in experiment. 

The remaining term of Eq. 10, S d , accounts for the 
additional spectroscopic data which can be included in the 
fit: lifetimes, branching ratios, multipole mixing ratios and 
known matrix elements. The summation extends over the 
number of such data points, rid, given for each data type, 
d , and user-defined weights, Wd, which are common for a 
given group of spectroscopic data. 

S d = Y J W d Y J -^( D n d - De n d )\ ( 13 ) 

d n d n d 

where D^ d and D^ d are the values of the spectroscopic 
data calculated using the current set of best-fit matrix 
elements and the experimental value, respectively. 

A simultaneous fit of a large number of free parameters 
(matrix elements), having unknown correlations and pos¬ 
sibly very different influences on the data, prevents pre¬ 
cise determination of degrees of freedom. In the simplest 
cases without these issues, the number of degrees of free¬ 
dom would be defined as a result of the subtraction of 
the number of experimental data points and the number 
of fitted parameters. The y 2 function resulting from the 
GOSIA calculations is normalised to the number of data 
points, including experimental intensities, branching ra¬ 
tios, lifetimes, mixing ratios and known matrix elements. 
In practical situations one deals exclusively with total y 2 
values, thus the normalised y 2 value yielding from the 
GOSIA code should be multiplied by the number of data 
points given, regardless of the user-defined weight, w. 

4.2.1 The GOSIA2 code 

When lifetimes of the lowest excited states are not known 
with sufficient precision, the measured Coulomb-excitation 
cross-sections need to be normalised in a different way, 
for example to the target excitation, as described in Sec¬ 
tion 4.1.3. The GOSIA2 code was developed to handle the 
simultaneous analysis of both target and projectile exci¬ 
tation. The x 2 function of Eq. 10 is minimised in parallel 


for the target and projectile whilst sharing the normali¬ 
sation factors as parameters across both functions. Using 
literature values of relevant matrix elements in the target 
nucleus, the normalisation constants can be constrained 
by the 7 -ray intensities of the target de-excitation. The 
solution then corresponds to the global minimum of the 
total x 2 function defined as the sum of x 2 functions for 
both reaction partners. If only two matrix elements are 
used to describe the excitation of the nucleus under study, 
a two-dimensional plot of the total x 2 surface may be used 
to evaluate uncertainties of fitted matrix elements, as de¬ 
scribed in more detail in Section 4.3.2. However, there 
are certain limitations of the code: when more unknown 
matrix elements are involved, estimation of their errors 
becomes more complicated and one of the procedures de¬ 
scribed in Sections 4.4.2 and 4.4.3 are required. 


4.3 Methods of error estimation 

4.3.1 Standard error estimation in GOSIA 

Statistical errors of the matrix elements are estimated af¬ 
ter the convergence of the global minimum of the x 2 func¬ 
tion and can be obtained from the probability distribution 
around the minimum. The applied method involves two 
steps. At first, the diagonal, or uncorrelated, uncertainties 
are calculated by sampling each matrix element about the 
minimum of the y 2 surface, finding the point where an in¬ 
crease in x 2 is achieved, satisfying the lcr condition. This 
condition is defined by requesting that the total integrated 
probability distribution in the space of the fitted param¬ 
eters be equal to the lcr confidence limit - 68.27 % [4]. 
At the same time, a multi-dimensional correlation matrix 
is built, which is then used in the second step in order 
to compute the fully correlated errors on each matrix ele¬ 
ment, satisfying the same condition. 


4.3.2 Two-dimensional y 2 surface analysis 

In a multi-parameter analysis, the global best fit can be 
found by constructing a y 2 hyper-surface with respect to 
all parameters. In the case of a two-parameter system one 
is able to visualise a 2 -dimensional y 2 surface as shown 
Figure 4(a). Here the example is of the two matrix el¬ 
ements, ( 2^11 JT 2 || 0 j'“) and ( 2 ^||U 2 || 2 j l ‘), usually sufficient 
to describe the excitation process of an even-even nucleus 
if only the 2+ —> Oj 1 " transition is observed. The minimum 
of such a surface, y 2 lin , can easily be found and the lcr- 
uncertainty contour can be defined as the region of the 
surface for which y 2 < y 2 lin + 1 [35-37]. This technique 
was used for the analysis of 94 ’ 96 Kr [38]. 

If one of the parameters is independently measured, 
e.g. via lifetime measurements, the y 2 surface can be eas¬ 
ily recalculated by adding the y 2 contribution of the new 
measurement at every point. A new lcr contour is then also 
defined, as shown in Figure 4(b). This goes too for other 
independent Coulomb excitation measurements, which 
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Fig. 4: (a) A full two-dimensional y 2 surface with respect 
to (2+11.52110+) and (2+||52||2+) for the 62 Fe projectile, 
reproduced from Ref. [39] . The data are normalised to the 
excitation of a 4.0-mg/cm 2 thick 109 Ag target at a beam 
energy of 2.86 MeV/u using GOSIA2. (b) The resulting 
surface when combined with lifetime measurements [40, 
41] and a ler cut applied with the condition that y 2 < 
Amin + 1- The individual lcr contours for the Coulomb- 
excitation and lifetime data are shown by the solid and 
dashed lines, respectively. Reproduced from Ref. [39] 


may come from the segmentation of a data set into an¬ 
gular ranges (see example in Figure 5) or different targets 
as described earlier. The final uncertainties are obtained 
by projecting the lcr uncertainty contour on the respective 
axes. While the projected uncertainties are useful for un¬ 
derstanding the precision on a given spectroscopic observ¬ 
able, such as B(E 2) values or spectroscopic quadrupole 
moments, the existing correlation between these parame¬ 
ters is lost. 

In the past, the assumption that the influence of the 
spectroscopic quadrupole moment, Q s , is negligible, or 
that otherwise its value can be assumed to be equal to 
zero has sometimes been used. In the case of 62 Fe, if it 
was not for the independent lifetime experiments, shown 
by the solid black line in Figure 4(b), the 2-dimensional 
lcr surface would not be constrained. One possibility is 
to project such a surface with a single value of Q s , or 
calculate only a 1 -dimensional surface at a fixed value of 
Q s . However, this would greatly underestimate the uncer¬ 
tainty, since correlations are ignored. In the 204 Rn example 
of Figure 5, this leads to a factor of 3.5 reduction in the 
true uncertainty. Instead, it would be preferable in these 



Fig. 5: A two-dimensional y 2 surface with respect to 
(2+||52||0+) and (2+||52||2+) for 202 Rn [42], A lcr cut is 
applied with the condition that x 2 < Xmin + 1 ■ The data is 
normalised to the excitation of a 4.0- mg/cm 2 thick 109 Ag 
target at a beam energy of 2.90 MeV/w using GOSIA2. The 
data was sub-divided into five different scattering angular 
ranges and their individual lcr limits are represented by 
the different bands; in increasing order of centre-of-mass 
scattering angle these are: solid black, dashed black, dot¬ 
ted black, solid grey and dashed gray. 


cases to use a model assumption where necessary to pro¬ 
vide limits of Q s as a function of (2+||52||0+), for example 
the rigid rotor model [43]. The total surface can then still 
be constrained but with a reasonable consideration of the 
uncertainty due to the influence of Q s . 

The graphical method however, becomes computation¬ 
ally time consuming and visually useless as the number 
of parameters increases. Therefore alternative solutions of 
error estimation are proposed and some examples are pre¬ 
sented in the following sections. Their applicability de¬ 
pends on the strength of the correlations between matrix 
elements. 


4.4 Selected applications 

4.4.1 Normalisation to the B(E 2) extracted from data sets 
where no correlations are observed 

The influence of the quadrupole moment of a given state 
on its excitation probability varies significantly with scat¬ 
tering angle as shown in Figure 1. This dependence can be 
exploited in order to determine both the transition prob¬ 
ability and the diagonal matrix element, even if only one 
7 -ray transition is observed in the nucleus of interest. If 
the particle detector covers a sufficiently broad range of 
centre-of-mass scattering angles, the simplest solution de¬ 
scribed in Ref. [44] can be applied. Here, in the first step, 
the B(E 2; 2+ —>■ 0+) value is derived from the excitation 
cross-section of the 2 ^ state for the lowest scattering an¬ 
gles. The influence of the quadrupole moment, Q s (2+), 
on the excitation probability of the 2 + state for this range 
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of scattering angles was estimated at 4% , which was be¬ 
low the statistical error of 7% of the corresponding 7 - 
ray intensity. It was therefore a reasonable approxima¬ 
tion to assume that in this case the observed transition 
strength depends only on the transitional matrix element. 
The adopted uncertainty of the B(E 2; 2+ —► 0^) included 
contributions from the statistical error of measured 7 -ray 
intensities in 44 Ar and 109 Ag, as well as the uncertainty 
on the relative 7 -ray efficiency, target matrix elements 
and the systematic error of 4% resulting from neglecting 
Q s (2+) in the Coulomb-excitation calculations. In the sec¬ 
ond step, this B(E2; 2^ —> Of) and its uncertainty were 
used in the further analysis as an additional data point 
in a GOSIA fit. The remaining data was then subdivided 
into three angular ranges, with the width and number of 
ranges being chosen to obtain the maximum sensitivity to 
Q S {2 J [). The 7 -ray intensities of 44 Ar from these ranges 
were normalised to the intensity measured for the first 
range, with relative normalisation factors fitted using the 
corresponding 109 Ag 7 -ray intensities. Then the standard 
version of the GOSIA code was used to simultaneously fit 
all the transitional and diagonal matrix elements to the 
measured intensities. 

4.4.2 Multiple Coulomb excitation and normalisation with a 
dominant transition to target excitation: combined 
GOSIA-GOSIA2 analysis 

In multiple Coulomb excitation of even-even nuclei, sev¬ 
eral states can be populated. In such cases the 2+ state 
is usually dominantly populated as compared to other 
excited states. When the lifetime of the 2^ state is not 
known with sufficient precision and the B(E2;2+ —>■ 0^) 
value cannot be extracted as described in Section 4.4.1, 
measured Coulomb-excitation cross sections need to be 
normalised in a different way using e.g., target excitation. 
However, a full analysis with the GOSIA2 code, as pre¬ 
sented in Section 4.4.3, is not possible as the number of 
parameters increases significantly. The error estimation in¬ 
cluding correlations between all matrix elements involved 
becomes very complex and practically impossible. A differ¬ 
ent solution needs to be found that handles both aspects: 
(i) normalisation to the target excitation and, (ii) error 
calculations including correlations between all matrix ele¬ 
ments. In such cases a combined analysis is required with 
the use of both standard GOSIA and GOSIA2 codes. 

In the first step, a simplified analysis is performed aim¬ 
ing to determine the B(E 2; 2+ —► Of) value for the projec¬ 
tile. Only one-step excitation of the 2f state is considered, 
taking into account that population of the 2 f state de¬ 
pends predominantly on both the B(E2;2f —► Of) value 
and spectroscopic quadrupole moment, Q s (2f). In order 
to gain sensitivity on the extraction of the quadrupole 
moment of the 2 f state, the data are divided in terms 
of particle-scattering angular range. The influence of the 
multi-step excitations resulting in population of higher- 
lying states is not usually included at this stage, although 
the level energies and a set of fixed “starting” matrix el¬ 
ements can be declared if reasonable assumptions can be 


made concerning their magnitudes and relative signs. The 
analysis is performed as described in Section 4.2.1 using 
the GOSIA2 code. As a result a two-dimensional total y 2 
surface (being the sum of y 2 for the projectile and target 
system) as a function of the B(E2\2f —► Of) value and 
the quadrupole moment Q s ( 2 f) is determined and reflects 
correlations between these two parameters. The final val¬ 
ues are determined by the minimum of the y 2 function 
and their error bars are obtained by projecting the la- 
contour on the respective axes, as in Section 4.3.2. The 
extracted B(E 2;2f —►Of) value is a first approximation. 
Its uncertainty includes: (i) the uncertainties of the 7 -ray 
intensities originating from the target excitation, (ii) the 
uncertainties of the 7 -ray intensities originating from the 
projectile excitation and, (iii) the uncertainties of the rel¬ 
evant, literature B(E 2) values for the target nucleus. 

In the second step, correlations with all remaining ma¬ 
trix elements, which couple higher-lying excited states ob¬ 
served in the experiment have to be investigated. This is 
performed using the standard GOSIA code with full error 
estimation procedure (see Section 4.3.1) implemented in 
GOSIA. All states populated in the Coulomb-excitation ex¬ 
periment, together with all observed 7 -ray intensities are 
taken into account in this part of analysis. All involved 
electromagnetic matrix elements are now introduced as 
well. Data extracted from the simplified GOSIA2 analysis, 
specifically (Of ||l£2||2f), serves as an absolute normali¬ 
sation for the standard GOSIA calculations. It is declared 
together with its uncertainty as an additional data point 
and thus it is treated in the fit on equal rights as the 7 -ray 
intensities. Other spectroscopic data i.e., 7 -ray branching 
ratios, mixing ratios, can also be included at this stage of 
analysis if known. Note that the (2f ||A2||2f) diagonal ma¬ 
trix element extracted for the projectile in the first part of 
the analysis is not included as an additional data point in 
the fit when switching to the standard GOSIA calculations. 
Information on (2f ||A2||2f) is implicitly given by the rel¬ 
ative normalisation constants extracted from the target 
excitation linking different angular data subdivisions. 

In order to link each data set resulting from subdi¬ 
vision into several particle-scattering angular ranges, the 
relative normalisation constants are required. These are 
usually calculated from the target excitation. The stan¬ 
dard GOSIA fit of observed 7 -ray intensities depopulating 
excited states in the target nucleus is performed using 
literature values of all relevant matrix elements (see Sec¬ 
tion 4.1.2). Calculated relative normalisation constants, 
Cij , for each data set are then further used to fit the 
projectile excitation. A small correction is applied here 
to achieve the same relative normalisation constants ob¬ 
tained in the GOSIA2 solution, where the projectile data 
is also considered. This is calculated using the ratio of 
the calculated yields for the normalisation transition in 
the target, from the GOSIA2 and standard GOSIA 

solutions. The fixed coupling of the relative normalisa¬ 
tion constants removes the knowledge of the uncertainty 
in their ratio. In order to preserve such information in the 
fit, it is included indirectly. The uncertainty of the 7 -ray 
intensities related to the 2f —► 0 + transition in the pro- 
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jectile, AI'P(2 J [ —> 0+), to which we are normalising, is 
defined in GOSIA so that it encompasses the uncertainty 
from the target excitation: 

. (M) 

where IH. and A'I? are the efficiency-corrected intensity 
and its associated uncertainty, of the 2+ —» 0+ transition 
in the projectile, respectively, and AI*. can be expressed 
as: 


AH, 


A%\ 2 (AB(E2-i ->■ g.s.)\ 

) V B{E2;i-> g.s.) ) 


(15) 


where 7* and A'l^. are the sum of efficiency-corrected 
intensities and associated uncertainties (in quadrature) of 
transitions depopulating a state i in the target, respec¬ 
tively. This assumes that this state is dominated by single- 
step excitation from the ground state and consequently by 
B(E2;i —► g.s.) and its uncertainty, AB(E2;i —► g.s.). 

As a result of the second part of the analysis with the 
use of the standard GOSIA code, a set of electromagnetic 
matrix elements between all states populated in the ex¬ 
periment is extracted. Note that the (0+||.E2||2+) matrix 
element, used as an absolute normalisation for the full 
standard GOSIA fit, originates from the simplified GOSIA2 
calculations where multiple Coulomb excitation was not 
necessarily correctly considered. This influence needs to 
be taken into account. For this purpose, the GOSIA2 cal¬ 
culations have to be repeated using the set of matrix ele¬ 
ments extracted in the second step of the analysis. Only 
(0+||i?2|12+) and (2+1)712112+) for the projectile are 
scanned as in the first approximation, while all the other 
matrix elements for the projectile are fixed and those for 
the target remain free. As a result, a new total \ 2 surface 
is calculated. Again, the (0+||712||2+) matrix element is 
determined from the \ 2 < Amin + 1 condition. It may dif¬ 
fer from the value obtained from the first approximation 
since the correlations with other matrix elements will be 
different. If this is the case, a full standard GOSIA analysis 
with the updated value of the (0+||i?2||2+) matrix ele¬ 
ment has to be repeated in order to achieve consistency. 
The whole standard GOSIA - GOSIA2 procedure should be 
iterated until the converged solution for both transitional 
and diagonal matrix elements for the 2+ state is obtained. 
A schematic procedure of the GOSIA - GOSIA2 analysis is 
presented in Figure 6. 

In some cases, such as 196 Po [26], the y-ray intensity 
of higher-lying transitions is too weak to be reliably ob¬ 
served in each of the angular subdivisions. For this data 
to be included, an additional data set must be declared 
in the GOSIA stage of the analysis that represents the sum 
total of all angular ranges. The simplest way to normalise 
this data is to use the total intensity of the 2+ —► Of 
normalisation transition allied with the B(E2) value that 
is already declared. This intensity then exclusively con¬ 
strains the absolute normalisation of the total data set 
with an uncertainty determined by the combination of the 



Fig. 6: A scheme of the combined analysis performed with 
the standard GOSIA and GOSIA2 codes. The presented 
method is used when normalisation to the target exci¬ 
tation is required in multiple-step Coulomb excitation of 
even-even nuclei. The red matrix elements in the level 
scheme of the figure are kept fixed during the GOSIA2 cal¬ 
culations, while the blue matrix elements are scanned to 
produce a 2-dimensional total y 2 surface plot (top right). 
All matrix elements are varied in the full GOSIA minimi¬ 
sation and the best fit values are used in the next GOSIA2 
calculation. Convergence is reached when the blue matrix 
elements are consistent in both GOSIA and GOSIA2 calcu¬ 
lations. 


B(E2) and 7 7 (2+ —> 0+) uncertainties. During the corre¬ 
lated error calculation, this uncertainty on the absolute 
normalisation is effectively propagated to the higher-lying 
transitions. 


4.4.3 Normalisation to target excitation when multiple 
single-step excitations are observed 

When multiple, single-step excitations are observed with 
similar intensity, such as in odd-mass systems, there are 
too many parameters to make an analysis of a full y 2 
hyper-surface feasible. Instead, a one-dimensional surface 
is constructed for each matrix element by scanning the pa¬ 
rameter to be investigated. At each point, the investigated 
parameter is kept fixed while all others are minimised 
with respect to \ 2 . For this, the minimisation procedure of 
GOSIA2, described in Section 4.2.1, is invoked. This proce¬ 
dure traces the lowest value path through the valley of the 
hyper-surface, effectively projecting the correlated surface 
to a given parameter. The constructed surface can then 
be used in order to extract the lcr uncertainty using the 
standard x 2 lin + 1 method [37]. There is an assumption 
here of parabolic behaviour about the minimum, which 
for strongly correlated systems may not necessarily be true 
and asymmetric limits may be obtained. 
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Computationally, the time involved to minimise the 
full parameter space hundreds of times is very large. For 
this reason, alternative methods of normalisation are pre¬ 
ferred, but Coulomb excitation of odd-mass or odd-odd 
systems with RIBs tend to lack the required lifetime and 
multipole mixing ratio data to sufficient precision. This 
approach has been successfully used for the analysis of 
Miniball experiments on odd-mass Sn isotopes [45] and 
the odd-odd 26 Na [46]. 

4.4.4 Normalisation to target excitation in a strongly 
correlated odd-mass system 

In the example of 9 'Rb, a strongly-deformed band built 
on the 3/2+ ground state is populated in Coulomb exci¬ 
tation with a 60 Ni target [47]. Mixed E2/M1 I —> I — 1 
transitions are roughly one order of magnitude stronger 
in intensity than I —> I — 2 transitions. In order to ex¬ 
tract transition probabilities in the low-energy part of the 
band, normalisation to target excitation is necessary. On 
the other hand, transition probabilities between the states 
that can only be reached in multi-step excitation are re¬ 
lated to measured intensity ratios in the nucleus of inter¬ 
est. As an example, the 4+ —»• 2+j2\ -A 0+ intensity ra¬ 
tio observed in Coulomb excitation of a weakly deformed 
even-even nucleus, assuming quadrupole moments equal 
to zero, depends exclusively on the (4] l '||i?2||2] 1 ") matrix 
element; changing the ( 0 +||.E 2 || 2 ^) matrix element would 
influence the total number of counts in both transitions, 
but not the ratio. This is no longer true if a significant frac¬ 
tion of nuclei (few percent) undergo excitation in each step 
of the process, which is often the case of deformed nuclei 
including 9 'Rb, but still observed relative intensities in the 
upper part of the band depend only weakly on lifetimes 
of the lowest-excited states. Therefore the analysis can 
be divided in two parts: the (7/2 + ||.E2||3/2 + ) matrix ele¬ 
ment is determined using normalisation to target excita¬ 
tion (gosia2 analysis), and the remaining matrix elements 
are extracted from the intensities measured for 9 'Rb us¬ 
ing the GOSIA code, fixing (7/2 + ||£’2||3/2 + ) at the value 
determined in the first part of the analysis. The choice of 
this matrix element was due to the fact that it corresponds 
to the only pure E2 transition from the ground state. The 
GOSIA2 code is used to find a minimum of the y 2 function 
resulting from comparison of measured and calculated 7 - 
ray intensities in 97 Rb and 60 Ni, as well as known spectro¬ 
scopic data in 60 Ni (B(E2;2 + —► 0 + ) and Q s (2+)). Mea¬ 
sured intensities of the 2 + —► 0 + transition in 60 Ni were 
scaled according to the measured beam composition, and 
their statistical uncertainties were adjusted to take into 
account the uncertainty of the beam composition, as de¬ 
scribed in Section 4.1.3. The minimisation is performed for 
several hundred starting values of (7/2 + ||£ l 2||3/2 + ) rang¬ 
ing from 0 to 3 eb. The (7/2 + ||i?2||3/2 + ) matrix element 
was fixed during the minimisation, while all other matrix 
elements are allowed to vary, with only constraints result¬ 
ing from Alaga rules [48] . In this way correlations between 
matrix elements are taken into account. Figure 7 presents 
the y 2 distribution as a function of (7/2 + ||7?2||3/2 + ) in 



Fig. 7: Total y 2 as a function of (7/2 + ||.E2||3/2 + ) in 9 'Rb. 
The open points show the y 2 obtained after convergence 
of the minimisation procedure, and the solid line is a poly¬ 
nomial fit of the y 2 distribution. The vertical lines corre¬ 
spond to y 2 nin (solid) and y 2 = y 2 lirl +1 (dashed; ler error 
bar). 


the vicinity of minimum. The vertical lines correspond to 
the adopted mean value (minimum of the y 2 distribution) 
and error bars (y 2 = xLin + 1) f° r the (7/2 + ||.E2||3/2 + ) 
matrix element. 

The second part of the analysis is performed using 
GOSIA with the (7/2 + ||£i2||3/2 + ) matrix element fixed at 
the value determined in the first part of the analysis. The 
errors of all remaining matrix elements are estimated us¬ 
ing the standard error evaluation procedure implemented 
in GOSIA (see Section 4.3.1). For transitions de-exciting 
states up to ll/ 2 + , it also is necessary to propagate the 
uncertainty of (7/2 + ||I?2||3/2 + ). For higher-lying transi¬ 
tions, contributions of this source of error to the total 
uncertainty is determined to be negligible. In this part 
of the analysis, (7/2 + ||i?2||3/2 + ) is fixed instead of be¬ 
ing fitted as an additional data point in order to make 
sure that its uncertainty is properly propagated. When 
(7/2 + ||.E2||3/2 + ), with the uncertainty determined in the 
first part of the analysis, is simply included in the fit on 
an equal basis as the 7 -ray transition intensities, its fi¬ 
nal uncertainty (and in consequence those of other ma¬ 
trix elements) is underestimated by the standard proce¬ 
dure of errors evaluation in GOSIA, as the y 2 minimum 
with respect to this matrix element is artificially made 
deeper by including two data points corresponding to the 
same observable ((7/2 + ||£’2||3/2 + ) and the 7/2+ —► 3/2+ 
transition intensities) in the fit. Such an effect has not 
been observed in the combined GOSIA-GOSIA2 analysis of 
Coulomb-excitation data in even-even nuclei (see 
Section 4.4.2) since there the 2+ — > g.s. transition is dom¬ 
inant, known with better precision than other 7 -ray inten¬ 
sities and thus its intensity serves basically to calculate the 
relative normalisation parameters for each experiment. In 
the case of 97 Rb, the 7/2+ -A 3/2+ transition is roughly 
20 times weaker than the strongest /—)•/— 1 transitions 
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Fig. 8: Total y 2 as a function of Q 0 in 97 Rb under assump¬ 
tion of the rigid rotor model for all E2 matrix elements. 
The open points show the \ 2 obtained after convergence of 
the minimisation procedure, the vertical lines correspond 
to Xmin (solid) and x 2 = Xmin + 1 (dashed; ler error bar). 
The weighted average of Q q values calculated from indi¬ 
vidual E2 matrix elements obtained in the full Coulex 
analysis including normalisation to the target excitation 
is shown in red. 


in this nucleus, which are consequently used to calculate 
the relative normalisation parameters. 


4.4.5 Normalisation to transition intensities in the nucleus 
of interest 

In very favourable cases of collective nuclei an estimation 
of transition probabilities can be obtained from the ratios 
of transition intensities in the nucleus of interest. It re¬ 
quires, however, strong model assumptions concerning the 
collectivity of the states (purely rotational or vibrational 
character). This procedure has been tested on the 97 Rb 
data where all E2 matrix elements between the observed 
states were coupled assuming the rigid rotor model. In 
this way one single parameter, corresponding to the tran¬ 
sitional quadrupole moment Q o of the band, was used to 
describe the E2 part of the measured gamma-ray inten¬ 
sities. No assumptions were made on the M1 matrix ele¬ 
ments of the mixed E2/M 1 / —> I— 1 transitions and in to¬ 
tal 7 parameters (one Q o value and 6 Ml matrix elements) 
were fitted to twenty measured y-ray intensities. In order 
to estimate the uncertainty of the extracted Qq value, the 
minimisation procedure was performed again for several 
hundred values of Q o kept fixed during minimisation with 
Ml matrix elements free to vary. A distinct minimum of 
the x 2 distribution was found as shown in Figure 8. Both 
the obtained value of Q o, as well as the error bars corre¬ 
sponding to x 2 = X'min+1 are consistent with the weighted 
average of Q o values calculated from individual E2 matrix 
elements obtained in the full Coulomb-excitation analysis 
including normalisation to the target excitation, presented 
in Sec. 4.4.4. 



Fig. 9: Example of a complicated detector shape in 9 and 
( t> coordinates: an off-centered annular detector with some 
parts damaged due to the high flux of incoming particles. 
The colors correspond to the number of detected events 
per pixel. 

4.5 Dealing with non-standard particle detectors 

Particle detectors used for RIB Coulomb-excitation ex¬ 
periments are usually axially symmetrical and have an 
absolute efficiency close to 100%. As long as the efficiency 
is uniform, any deviations from 100% are included in the 
normalisation constants (see Section 4.1). However, with 
the expected increase of RIB intensities, the standard an¬ 
nular Si detectors that are currently used will likely be 
replaced by more complicated particle detection set-ups, 
possibly consisting of various types of detectors differing 
in efficiency. In addition, radiation damage may deterio¬ 
rate parts of a detector, resulting in a very complicated 
shape in the 9-cf) plane. 


4.5.1 Complex particle-detector shapes 

In the example of 44 Ar [44], the beam was not well fo¬ 
cused and had a halo of about 0.5% of the total inten¬ 
sity, hitting the particle detector directly. Some parts of 
the particle detector had to be excluded from the analysis 
due to deterioration caused by the direct beam and re¬ 
sulting impossibility of distinguishing between direct and 
scattered beam. Together with a non-axial position of the 
beam spot, this resulted in a complicated shape of the de¬ 
tector in the 9-(f) plane (see Figure 9), which had to be 
taken into account during the Coulomb-excitation analy¬ 
sis using the standard GOSIA code. The standard methods 
of describing the particle detection geometry provided by 
the code did not allow a proper handling of this complica¬ 
tion. Therefore a new method was introduced and tested. 
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Each of the 1536 pixels of the detector (96 strips by 16 
annular rings) was approximated by a small circular de¬ 
tector, whose size was chosen to optimally reproduce both 
the absolute Rutherford cross-section and the calculated 
correction factors 1 for both 44 Ar and 109 Ag. The results of 
such approximation as compared to a standard integration 
procedure were verified for each strip and the differences 
were below 2% for all transitions. The correction factors 
depend weakly on actual values of matrix elements and 
thus the verification performed for the initial set of matrix 
elements remain valid throughout the minimisation pro¬ 
cedure. The difference between the quadrupole moment 
of the state obtained from the analysis with a proper 
detector shape taken into account and of that when it 
was assumed to be axially symmetric with all segments 
working, was around 20 %. 


4.5.2 Non-uniform particle-detector efficiency 

If the efficiency of the particle detector changes as a func¬ 
tion of scattering angle, this information should be in¬ 
cluded in the detector description used by the GOSIA code. 
This can be done by modifying the shape of the particle 
detector with respect to its real angular coverage. The 
simplest solution, used in the analysis of Coulomb excita¬ 
tion of 152 Sm [49] is to reduce the detector coverage in the 
(f> plane according to its relative efficiency. The 136 Xe ions 
scattered on the 152 Sm target were identified in the focal 
plane of the VAMOS spectrometer placed at 35°, which 
corresponds to the detection of ions scattered at angles be¬ 
tween 28° and 42° in 9 and -7° and 7° in </>. The simulated 
detection efficiency as a function of 9 scattering angle [50] 
is presented in 10 (a), and resulting particle detector shape 
included in GOSIA in 10(b): in the maximum of the effi¬ 
ciency curve the real coverage in (j> has been assumed, and 
for other scattering angles it has been scaled according to 
the efficiency. 

Such a solution works well if the effects of particle- 
7 -ray correlations can be neglected, i.e. when the 7 -ray 
detection set-up consists of many detectors placed sym¬ 
metrically in 9 and <f> and the 7 -ray intensities from all de¬ 
tectors are summed together. The efficiency curve should 
also be relatively smooth and simple, which is the case 
of the presented example. In other cases, however, such a 
modification of the particle detector shape may affect the 
calculated particle-gamma angular distributions and, in 
consequence, the extracted matrix elements. An alterna¬ 
tive method has therefore been tested, similar to the one 
presented in Section 4.5.1. The detector has been approx¬ 
imated by a set of 729 small circular particle detectors. In 

1 The correction factors, introduced in the GOSIA code in 
order to speed up the minimisation process, are defined for each 
observed 7 -ray transition as a ratio of its intensities calculated 
for a given set of matrix elements: the one integrated over the 
angular range covered by the particle detector and the range of 
incident energies resulting from slowing down of beam particles 
in the target, and that calculated for mean values of beam 
energy and scattering angle specified by the user. 



28 30 32 34 36 38 40 42 44 


6 [deg.] 

Fig. 10: Two methods to take into account non-uniform 
particle detector efficiency in GOSIA analysis, (a) Absolute 
efficiency of the particle detector as a function of 9 scatter¬ 
ing angle [50]. (b) Detector shape resulting from relating 
its coverage in </> to the efficiency; dashed lines correspond 
to the true coverage of the detector, (c) Approximation of 
the detector by a large set of pixel-like circular detectors, 
which sizes reflect the efficiency. 


the first step the particle detector was assumed to have a 
uniform 100 % efficiency, which corresponded to a rectan¬ 
gle in the 9-cj) plane or alternatively to all pixels having the 
same size. This size was adjusted to reproduce both the 
Rutherford cross section and correction factors for 152 Sm 
calculated for the rectangular particle detector. In the sec¬ 
ond step the size of each pixel was scaled according to the 
relative efficiency, as presented in 10 (c). 

The results of the two approaches were compared and 
were compatible within 2 % for excitation of states up to 
12 + . On the other hand, when the reduction of efficiency 
at the edges of the detector was neglected, the calculated 
relative 7 -ray yields differed by up to 14% as compared 
to that calculated taking the non-uniform efficiency into 
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account. The effect was the strongest for multi-step exci¬ 
tation and non-yrast states. 

4.6 Other sources of systematic errors 

Numerous approximations are used in the GOSIA code, de¬ 
scribed in details in Ref. [4]. They may amount to up to 
5% of the calculated 7 -ray intensity and thus very small 
error bars that may result from GOSIA error estimation 
procedure should be treated as statistical errors only and 
further adjusted to take into account the systematic er¬ 
rors. 

The most important source of systematic error is usu¬ 
ally related to the semiclassical approximation of the 
Coulomb-excitation process used in the analysis. This sim¬ 
plified treatment is expected to differ from a full Coulomb- 
excitation calculation by a factor of 1 /iv, where v is the 
Sonnnerfeld parameter [51], which for heavy ions (v ~ 10 2 ) 
amounts to a few percent. Other sources of systematic er¬ 
rors arising from approximations used in the GOSIA code 
are discussed for example in Ref. [13, 52] and most of 
them (corrections due to atomic screening, vacuum polar¬ 
isation, relativistic effects, E4 excitation) are found to be 
negligible. 

The deorientation effect (modification of the nuclear 
state alignment due to the interaction with the rapidly 
fluctuating hyperfine fields of the deexciting atom recoil¬ 
ing in vacuum) influences the 7 -ray angular distributions 
observed in Coulomb-excitation experiments. Current 
atomic model predictions of the deorientation effect are 
too computer intensive, and not yet viable, to consider 
their incorporation into GOSIA. Instead a schematic two- 
state model has been adopted with parameters fitted to 
available deorientation effect data [4, 53, 54]. Extensive 
studies [55] of the efficacy of the deorientation correction 
implemented in GOSIA have shown that, on average, the 
default values adopted in GOSIA work surprisingly well. 
On the other hand, averaging over particle and 7 -ray de¬ 
tection angles washes out sensitivity to the angular cor¬ 
relation effects for 7 -ray deexcitation and thus minimises 
the influence of deorientation on the results. In the in the 
104 Ru case [52] , changing the parameters of the deorienta¬ 
tion model by 20 % produced less than 2 % change in fitted 
matrix elements. 

The effect of virtual excitation of the giant dipole reso¬ 
nance can influence the excitation of low-lying states. This 
is taken into account using the concept of dipole polariz¬ 
ability [51] and applying a correction to the quadrupole 
interaction. This effect becomes important for light nuclei. 
In the analysis of 10 Be [56] it was found that the uncer¬ 
tainty of 25% on the polarizability parameter translated 
into 20 % uncertainty on the diagonal matrix element of 
the first-excited state. 

The integration procedures used in GOSIA to account 
for beam stopping in the target and the angular coverage 
of the particle detector may be quite sensitive to user- 
defined meshpoints. This is true in particular for compli¬ 
cated shapes of the particle detector, large ranges of in¬ 
cident energies (i.e. ’’thick-target” measurements, where 


beam is stopped in the target), small scattering angles 
and high energies of excited states (above 1 MeV in a 
single step). The influence of meshpoints on calculated in¬ 
tegrated 7 -ray intensities should be in any case verified 
and, if the differences between the calculated integrated 
yields for different sets of meshpoints are comparable with 
statistical uncertainties of the 7 -ray yields, should be in¬ 
corporated in these. 

Especially for well-deformed, or on the contrary, non¬ 
collective nuclei the lifetimes of Coulomb-excited states 
may be as long as nanoseconds. In such cases it is essential 
to take into account the modification of 7 -ray efficiency 
due to the modified solid angle covered by the 7 -ray de¬ 
tectors. This effect was observed for example in analyses 
of 9 'Rb [47] and 98 Sr [57] MINIBALL data and the af¬ 
fected transition intensities were either excluded from the 
analysis [57] or their statistical errors increased to take 
into account the modified efficiency [47]. 

The standard minimisation procedure works best if 
only E 2 matrix elements are needed to describe the ob¬ 
served excitation. The probability of getting trapped in a 
local minimum increases with every multipolarity included 
in the calculations. In particular, it is often observed that 
the errors on Ml matrix elements are underestimated. 
Many sets of starting values of matrix elements, includ¬ 
ing relative signs, should be tested before final values of 
matrix elements and their uncertainties are determined. 


5 Summary and outlook 

In summary, we have presented a number of methods 
for normalisation of Coulomb excitation data with Ra¬ 
dioactive Ion Beams (RIBs), using the GOSIA and GOSIA2 
codes. Analysis techniques have been presented with ref¬ 
erence to specific cases where the techniques were pio¬ 
neered. While excited-state lifetimes, in combination with 
other independent spectroscopic data, provide the sim¬ 
plest method of normalising Coulomb-excitation data, we 
have shown that it is possible to treat data in differ¬ 
ent ways, such as normalising to target excitation. These 
methods and techniques will gain an even greater im¬ 
portance as a wider range of post-accelerated RIBs be¬ 
come available at the next generation of ISOL facilities, 
such as HIE-ISOLDE [58], SPIRAL2 [59], ARIEL [60], 
CARIBU [61] and SPES [62]. In particular, the higher 
beam energies offered for heavy exotic nuclei will produce 
data for which multiple-step Coulomb excitation of iso¬ 
topes with a lack of spectroscopic data in the literature, 
becomes standard fare. 
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